Low-elevation glaciers in coastal regions of Alaska, the Canadian Arctic, individual ice caps around the Greenland ice sheet, and the Patagonia Ice Fields have an aggregate glacier area of about 332 × 103 km2 and account for approximately 42% of all the glacier area outside the Greenland and Antarctic ice sheets. They have shown volume loss, especially since the end of the 1980s, increasing from about 45% in the 1960s to nearly 67% in 2003 of the total wastage from all glaciers on Earth outside those two largest ice sheets. Thus, a disproportionally large contribution of coastal glacier ablation to sea level rise is evident. We examine cumulative standardized departures (1961–2000 reference period) of glacier mass balances and air temperature data in these four coastal regions. Analyses indicate a strong association between increases in glacier volume losses and summer air temperature at regional and global scales. Increases in glacier volume losses in the coastal regions also coincide with an accelerated rate of ice discharge from outlet glaciers draining the Greenland and West Antarctic ice sheets. These processes imply further increases in sea level rise.
Introduction
Acceleration of wastage of mountain and subpolar glaciers has been reported for many regions of the globe (Arendt et al., 2002; Kaser and Osmaston, 2002; Rignot et al., 2003; Meier et al., 2003). Glaciers in coastal regions have experienced exceptionally large volume losses, and glacier retreat during the late 1980s and early 1990s has been observed in Alaska, the Patagonia Ice Fields (PIF), coastal sections of large islands in the Canadian Archipelago, and glaciers around the Greenland Ice Sheet (Rignot et al., 2003; Thomas et al., 2003; Bradley and Serreze, 1987; Braun et al., 2001). We have found that amplification of glacier wastage appears to be related primarily to climate warming (Fig. 1). The purpose of this paper is to provide new quantitative evidence on the association between increases in summer temperatures and increases in volume losses of glaciers in several coastal regions.
Data and Methods
In this paper we analyze new and recently updated observations of glacier mass balance for a number of glaciers across the globe for the period 1960 through 2003 (Dyurgerov, 2002, updated in 2005). The data are from a number of sources (IASH (ICSI)–UNESCO, 1973, 1977, 1985; IAHS (ICSI)–UNEP–UNESCO, 1988, 1993, 1998, 2005; Cogley, 2002; Dyurgerov, 2002; Dyurgerov and Meier, 2005). We focus on glaciers in coastal regions located in four regions: Alaska (glacier area is ∼90 × 103 km2), coastal sections of large Islands in the Canadian Archipelago (glacier area is ∼151.8 × 103 km2), glaciers around the Greenland Ice Sheet (glacier area is ∼70 × 103 km2), and the Patagonia Ice Fields (glacier area is ∼19.9 × 103 km2). Different, but commonly accepted methods have been used to measure, calculate, and estimate accuracy in glacier mass balance changes in these regions (Mayo et al., 1972; Østrem and Brugman, 1991; Arendt et al., 2002; Dyurgerov, 2002).
To compute a time series of global annual mass balance (b30) we used mass balance time series for all glaciers with record lengths of 30 yr and longer (Fig. 1). For Alaska (Fig. 2B), we used observations from three benchmark glaciers (i.e., Gulkana, Wolverine, and Lemon Creek). The average annual mass balance time series for these three glaciers matches well the time series derived using repeated laser altimetry methods (Arendt et al., 2002; Meier and Dyurgerov, 2002). Data for the Canadian Archipelago were obtained from a number of sources (Cogley, 2002; Dyurgerov, 2002; data from R. Koerner, private communication).
We used all available observational mass balance data for Greenland [IAHS (ICSI)–UNESCO, 1985; IAHS (ICSI)–UNEP–UNESCO, 1988; Hasholt, 1988; Knudsen and Hasholt, 2003]; however, there are limited and sporadically measured data for these sites. These limited observations have been compiled during the last three decades [IAHS (ICSI)–UNEP–UNESCO, 1988; Hasholt, 1988; Dyurgerov, 2002; Knudsen and Hasholt, 2003] and indicate significant correlation with mass balances in western Svalbard and the Canadian Archipelago (correlation coefficient of 0.72, significant at a 0.95 confidence level). Because of these significant correlations, observations for western Svalbard and the Canadian Archipelago were used to reconstruct mass balances for the Greenland glaciers back to 1961 (Fig. 2C).
Because glacier mass balance data for the Patagonia Ice Fields (Fig. 2D) are extremely scarce, we used measurements of glacier thickness from surveys made in 1995 and 2000 by the Shuttle Radar Topography Mission (SRTM) (Rignot et al., 2003). These measurements have been compared with digital-elevation-model–derived topographic maps developed in 1975, and surface changes for 63 glaciers have been determined (Rignot et al., 2003).
Monthly temperature reanalysis data set from the National Center for Environmental Prediction (NCEP)/National Center for Atmospheric Research (NCAR) ( http://www.cdc.noaa.gov/cdc/reanalysis/reanalysis.shtml) was used to compute global and regional time series of annual and summer (June through August for the Northern Hemisphere, and December through February for the Southern Hemisphere) temperatures.
All time series (except data for the Patagonia Ice Fields) were standardized for ease of comparison using data for 1961–2000 as the reference period for the computation of departures. The standardization was performed by subtracting the long-term (i.e., 1961–2000) mean of a time series from each value and then dividing by the respective long-term standard deviation. Subsequently, the time series of standardized departures were integrated to produce time series of cumulative standardized departures. Data for the Patagonia Ice Fields were not standardized because only three estimates of annual glacier mass balance were available. These data were included in the analysis because they provide some indication of glacier mass changes for glaciers in the Southern Hemisphere.
Global and regional time series of cumulative standardized departures of annual glacier mass balance and summer air temperature were compared to identify associations between these time series and to identify coincident shifts in the time series.
Acceleration in Glacier Wastage in Coastal Regions
Comparison of annual glacier mass balance with summer temperature for the Canadian Archipelago indicates a close association between these time series. In addition, an increase in the wastage of Canadian Archipelago glaciers during the late 1980s appears to coincide with an increase in summer temperature (Fig. 2A). Increases in summer temperature and melting have resulted in upward movements of the equilibrium line altitude (ELA) and runoff line altitude. Ice caps in the Canadian Archipelago are especially sensitive to changes in the ELA due to their geometry (Bradley and Serreze, 1987; Braun et al., 2001; Burges and Sharp, 2004). Similarly, time series of annual glacier mass balance for Alaska indicate an increase in mass wastage during the late 1980s that also appears related to an increase in summer air temperature (Fig. 2B). This result is supported by changes in mass balance and glacier volume in Alaska derived from laser altimetry surveys (Arendt et al., 2002). These surveys indicate unprecedented acceleration of glacier wastage for glaciers in Alaska since the late 1980s (Fig. 2B).
The reconstructed mass balances for Greenland also show a shift toward increased mass wastage. However, the acceleration in mass wastage for Greenland appears to have occurred later than the shift observed for the Canadian Archipelago and Alaska, i.e., around 1997, coincident with increases in summer air temperature during that time. An acceleration of glacier volume losses in Greenland is consistent with rapid surface thinning of ice observed in the southwestern and southeastern parts of the Greenland Ice Sheet; this thinning has propagated to elevations as high as 1500 and 2000 m a.s.l. (Thomas et al., 2003; Zwally et al., 2002; Krabill et al., 2004).
Other evidence of increased glacier mass wastage is the increased loss of mass from outlet and tidewater glaciers in the Northern and Southern Patagonia Ice Fields. Figure 2D suggests that glacier mass balance has declined rapidly in the Patagonia Ice Fields during the past few decades. It has been estimated that Patagonia Ice Fields have experienced a volume loss of 41.9 ± 4.4 km3 yr−1 during 1975–2000 (Rignot et al., 2003). In addition, the average rate of thinning of this ice has more than doubled during the 1995–2000 period (Rignot et al., 2003). There is growing evidence that glacier retreat and volume loss in the Patagonia Ice Fields also has accelerated (Fig. 2D), but the cause of the acceleration remains unknown because of a scarcity of meteorological records on or near the Patagonia Ice Fields. There may be several causes of the apparent accelerated glacier wastage, such as increase in thinning rate in ablation areas, resulting from longitudinal stretching, ice calving into lakes, and other dynamic processes related to increased melt water on the surface and on the glacier bed (Rivera and Casassa, 2004). Warming during the winter and summer seasons during about the past 40 yr has been approximately 0.5°C at the 850 hecto-Pascal atmospheric pressure level. This warming has resulted in about a 5% decrease in the percentage of precipitation that occurs as snow and also has increased the annual melt rate in ablation areas by about 50 cm yr−1 (Rasmussen et al., 2003; Carasco and Quintana, 2003). In addition, other researchers have reported a slight increase in winter precipitation and increases in both winter and summer temperatures during recent years for the southern part of South America and the Antarctic Peninsula (Carasco and Quintana, 2003). The winter warming was found to be larger than the summer warming in these regions. Examination of temperature data from the NCEP/NCAR reanalysis data set indicates increased winter temperatures near the Patagonia Ice Fields after the late 1980s, but there is no increase evident in summer temperatures. Because of few observations of temperature near the Patagonia Ice Fields and the low resolution of the reanalysis data, these conclusions are only suggestive.
An important and common feature of near-coastal glaciers is that they are at low elevations, and their frontal parts, grounded or floating, are near sea level where increases in air temperature may have an enhanced effect on the melting of snow and ice. Previously it has been shown that the magnitude of summer melting is non-linearly related with increases in summer air temperature. A cubic parabola is a reliable approximation for the relation between glacier-summer melt and June–August air temperature (Krenke and Khodakov, 1966). Thus, low-elevation glaciers are highly sensitive to increases in air temperature (Meier, 1984; Oerlemans, 2001). In addition, the ratio of rain to snow increases with increases in air temperature, and increases in this ratio may have strong effects on melting because of reductions in average surface albedo and increases in absorbed solar radiation. Although glaciers in near-coastal regions generally occur in places where precipitation is high, these glaciers have been found to be mostly sensitive to changes in air temperature (Naito et al., 2001).
The amplification of melting rates for near-coastal regions has affected glaciers in coastal sections of southwestern and northwestern Scandinavia, as well. This is somewhat unusual in that these glaciers have accumulated mass since the start of observations in the early 1960s and during recent decades, when global temperatures have increased and glacier wastage has accelerated for most other glaciers worldwide (Meier et al., 2003). However, glaciers in the coastal areas of Scandinavia have started to lose mass and retreat during the past few years, possibly responding to observed changes in temperature. The recent observations of glacier mass changes for this region suggest that changes in Scandinavian glaciers are now consistent with the changes in mass balances observed for other coastal and continental regions around the world during recent decades (Meier et al., 2003; Kjøllmoen, 2003a, 2003b, 2004).
Oceanic Forcing
Another major forcing of mass changes for near-coastal glaciers is water temperature. Increases in global sea-surface temperature (0–3000 m depth) during 1955–2003 have been reported recently (Levitus et al., 2005). At the base of glaciers where floating ice is in contact with oceanic water, glacier melting may have exceeded the rate of surface melting by several orders of magnitude. Bottom melting along the 50 km floating length of the Petermann outlet glacier in Greenland (81°N) has been estimated to be on average 7 m yr−1, reaching in some places as much as 50 m yr−1 (Steffen et al., 2004). Similarly, melting along the coastal sections of the Greenland ice sheet has been estimated to be approximately 10 m yr−1 (Thomas, 2004), and up to 5.5 m yr−1 for the West Antarctic ice sheet (Shepherd et al., 2004). Observations of mass loss for coastal regions of Greenland and West Antarctica suggest that the changes in mass and the driving processes are similar to the processes occurring in regions outside the ice sheets. Jakobshavn outlet glacier, which drains about 7% of the entire Greenland Ice Sheet, indicates an acceleration of flow rates and mass wastage since about the early 1990s (Thomas et al., 2003; Joughin et al., 2005). In addition, rapid disintegration of ice shelves and tidewater glaciers in West Antarctica and in the High Arctic has been observed (Scambos et al., 2000; Mueller et al., 2003). Ward Hunt Ice Shelf, which was the largest ice shelf in the Arctic, has rapidly disintegrated during recent summers, most likely due to increases in surface melting and the warming of sea water (Mueller et al., 2003).
Glacier Sensitivity to Temperature
A comparison of changes in glacier mass balance with changes in summer temperature (Table 1) shows the sensitivity db/dTs of glacier mass balances to changes in summer air temperature (here b is the annual mass balance and Ts is the summer air temperature averaged for three summer months). For the period 1961–2001, our observed globally averaged db/dTs was −380 mm °C yr−1, which is reasonably close to the −350, −370, −390, and −410 mm °C yr−1 calculated by different authors and methods (see Raper and Braithwaite, 2006). These global-scale similarities imply that serious progress has been made in glacier mass balance observations and modeling.
We have recently shown that the sensitivity of glaciers to temperature in the Northern Hemisphere has both increased and become dramatically less variable since 1987 (Dyurgerov, 2006). We suggest that these changes in sensitivity and its variability may have resulted from changes in zonal Ts distribution patterns (Fig. 3). Before the late 1980s temperature anomalies across the Earth were both positive and negative (Fig. 3). Since the late 1980s zonally averaged annual air temperature anomalies have been positive for almost all latitudinal zones (Figs. 1 and 3). These changes are global phenomena, likely related to global climate change but not yet fully explained. Glaciological observations show that these changes affected glaciers more at low elevation and near the ocean. The acceleration of volume losses for glaciers in near-coastal regions is likely indicative of larger changes at the margins of the Greenland and Antarctic Ice Sheets. Recent observations, during the 1991 through 2000 period, indicate large volume losses of the Greenland Ice Sheet that have averaged approximately 78 km3 yr−1. These volume losses correspond to a 0.21 mm rise in sea level (Box et al., 2004). In addition, glacier thinning rates in some coastal regions of West Antarctica (e.g., Amundsen Sea) for recent years are double those observed during the 1990s (Thomas et al., 2004). Vaughan (2006) suggests that mass loss from the Antarctic Peninsula as a direct response to global warming has been significant (i.e., 0.008–0.055 mm yr−1) and could triple over the next 50 yr.
Increase in Glacier Wastage in Coastal Regions and Sea Level
Several recent studies have provided new estimates of sea-level rise (SLR) that are not consistent with previously published values (e.g., Church and Gregory, 2001). Cabanes et al. (2001) suggest that previously reported rates of SLR have been 2–3 times too high because many of the gauges that measure SLR are located in coastal areas where ocean temperatures are relatively warm. Cabanes et al. suggest that the true value for the thermosteric component is about 0.5–1.0 mm yr−1 (for 2001). Lombard et al. (2004) report an even smaller value of 0.3 mm yr−1 for the steric component and report a much larger residual value of SLR (observed minus thermosteric) of 1.4 ± 0.6 mm yr−1 (Lombard et al., 2004). Antonov et al. (2002) report a similar value of 1.4 mm yr−1, which was derived based on estimates of the amount of freshwater needed to increase sea level. Antonov et al. suggest that the ice-melt on land, mostly from glaciers, is the source of this freshwater. After a reanalysis of tide-gauge data, Miller and Douglas (2004) stated that evidence suggests that changes in ocean temperature and salinity only account for a fraction of SLR and that the change in mass (the eustatic component) played a dominant role in 20th century global SLR. The findings of Miller and Douglas are consistent with the non-thermosteric value of SLR of 1.4 mm yr−1 by Antonov et al. (2002) and Lombard et al. (2004).
The increased wastage observed for glaciers in near-coastal regions has resulted in an increased contribution of near-coastal glaciers to globally averaged glacier volume loss (Fig. 4). The acceleration of wastage observed for near-coastal glaciers also is likely in progress at the largest ice sheets, which suggests that larger than previously estimated eustatic contributions to SLR may be occurring.
Conclusion
We compiled new mass balance observations for glaciers in coastal regions, filled gaps in some time series, and extended time series for Alaska, the Canadian Archipelago, and individual ice caps around the Greenland Ice Sheet back to 1961. We also synthesized new and previously published glacier mass balance observations. Analysis of these records indicates an increase in ice wastage for near-coastal glaciers since the late 1980s or early 1990s. This volume loss has increased from about 45% in the 1960s to 67% in 2003 of the total mass wastage from all glaciers on Earth outside the two largest ice sheets. Thus, a disproportionally large contribution of coastal glaciers to sea level is evident. These increases in volume losses are strongly associated with increases in surface summer air temperature. In addition, near-coastal glaciers indicate an increased sensitivity to air temperature. The increase in temperature sensitivity appears to have occurred during 1987–2001 for glaciers in the Canadian Archipelago, around 1994–2001 for glaciers in Alaska, during 1998–2001 for Greenland ice caps, and near 1996–2001 for the Patagonia Ice Fields. The large negative mass balances for glaciers in coastal regions are in agreement with reports of an acceleration in surface and basal melting, iceberg calving in Greenland, and observed disintegration of ice shelves in West Antarctica.
A concern raised by these results is that the acceleration in wastage of mountain glaciers and subpolar ice caps in coastal areas may be an early warning of larger changes in two ice sheets; Greenland and West Antarctic. In addition, recent observations suggest that changes of these large ice sheets may have already started.
Acknowledgments
We thank Mark Meier and John Hollin for many important comments and suggestions that helped to improve the paper. This study was supported, in part, through NASA grant NNG04GM09G for the GLIMS Core Functions (R. L. Armstrong, Principal Investigator), and by NASA grant NAG5-13691 and U.S. National Science Foundation grant NSF/OPP-0425488 (M. Dyurgerov, Principal Investigator).